Statistical analysis of neuronal data: Development of quantitative frameworks and application to microelectrode array analysis and cell type classification

With increasing amounts of data being collected in various fields of neuroscience, there is a growing need for robust techniques for the analysis of this information. This thesis focuses on the evaluation and development of quantitative frameworks for the analysis and classification of neuronal data from a variety of contexts. Firstly, I investigate methods for analysing spontaneous neuronal network activity recorded on microelectrode arrays (MEAs). I perform an unbiased evaluation of the existing techniques for detecting ‘bursts’ of neuronal activity in these types of recordings, and provide recommendations for the robust analysis of bursting activity in a range of contexts using both existing and adapted burst detection methods. These techniques are then used to analyse bursting activity in novel recordings of human induced pluripotent stem cell-derived neuronal networks. Results from this review of burst analysis methods are then used to inform the development of a framework for characterising the activity of neuronal networks recorded on MEAs, using properties of bursting as well as other common features of spontaneous activity. Using this framework, I examine the ontogeny of spontaneous network activity in in vitro neuronal networks from various brain regions, recorded on both single and multi-well MEAs. I also develop a framework for classifying these recordings according to their network type, based on quantitative features of their activity patterns. Next, I take a multi-view approach to classifying neuronal cell types using both the morphological and electrophysiological features of cells. I show that a number of multi-view clustering algorithms can more reliably differentiate between neuronal cell types in two existing data sets, compared to single-view clustering techniques applied to either the morphological or electrophysiological ‘view’ of the data, or a concatenation of the two views. To close, I examine the properties of the cell types identified by these methods.Supported by a Wellcome Trust PhD Studentship and a National Institute for Health Research (NIHR) Cambridge Biomedical Research Centre Studentshi

Comparison of burst detectors for spike trains

Accurate identification of bursting activity is an essential element in the characterization of neuronal network activity. Despite this, no one technique for identifying bursts in spike trains has been widely adopted. Instead, many methods have been developed for the analysis of bursting activity, often on an ad hoc basis. Here we provide an unbiased assessment of the effectiveness of eight of these methods at detecting bursts in a range of spike trains. We suggest a list of features that an ideal burst detection technique should possess and use synthetic data to assess each method in regard to these properties. We further employ each of the methods to reanalyze microelectrode array (MEA) recordings from mouse retinal ganglion cells and examine their coherence with bursts detected by a human observer. We show that several common burst detection techniques perform poorly at analyzing spike trains with a variety of properties. We identify four promising burst detection techniques, which are then applied to MEA recordings of networks of human induced pluripotent stem cell-derived neurons and used to describe the ontogeny of bursting activity in these networks over several months of development. We conclude that no current method can provide "perfect" burst detection results across a range of spike trains; however, two burst detection techniques, the MaxInterval and logISI methods, outperform compared with others. We provide recommendations for the robust analysis of bursting activity in experimental recordings using current techniques.Experimental data collection was supported by the BBSRC (PC, OP, grant number BB/H008608/1). EC was supported by a Wellcome Trust PhD Studentship and NIHR Cambridge Biomedical Research Centre Studentship. CWT was supported by a bursary from the Bridgwater Summer Undergraduate Research programme.This is the final version of the article. It first appeared from the American Physiological Society via https://doi.org/10.1152/jn.00093.201

Quantitative differences in developmental profiles of spontaneous activity in cortical and hippocampal cultures.

BACKGROUND: Neural circuits can spontaneously generate complex spatiotemporal firing patterns during development. This spontaneous activity is thought to help guide development of the nervous system. In this study, we had two aims. First, to characterise the changes in spontaneous activity in cultures of developing networks of either hippocampal or cortical neurons dissociated from mouse. Second, to assess whether there are any functional differences in the patterns of activity in hippocampal and cortical networks. RESULTS: We used multielectrode arrays to record the development of spontaneous activity in cultured networks of either hippocampal or cortical neurons every 2 or 3 days for the first month after plating. Within a few days of culturing, networks exhibited spontaneous activity. This activity strengthened and then stabilised typically around 21 days in vitro. We quantified the activity patterns in hippocampal and cortical networks using 11 features. Three out of 11 features showed striking differences in activity between hippocampal and cortical networks: (1) interburst intervals are less variable in spike trains from hippocampal cultures; (2) hippocampal networks have higher correlations and (3) hippocampal networks generate more robust theta-bursting patterns. Machine-learning techniques confirmed that these differences in patterning are sufficient to classify recordings reliably at any given age as either hippocampal or cortical networks. CONCLUSIONS: Although cultured networks of hippocampal and cortical networks both generate spontaneous activity that changes over time, at any given time we can reliably detect differences in the activity patterns. We anticipate that this quantitative framework could have applications in many areas, including neurotoxicity testing and for characterising the phenotype of different mutant mice. All code and data relating to this report are freely available for others to use.PC and AM were supported by the Wellcome Trust Genes to Cognition programme. PC received additional support from the Biotechnology and Biological Sciences Research Council (BB/H008608/1). EC was supported by a Wellcome Trust PhD studentship and Cambridge Biomedical Research Centre studentship. SJE was supported by an Engineering and Physical Sciences Research Council grant (EP/E002331/1).This is the final published version. It first appeared at http://link.springer.com/article/10.1186%2Fs13064-014-0028-0

Characterization of Early Cortical Neural Network Development in Multiwell Microelectrode Array Plates.

We examined neural network ontogeny using microelectrode array (MEA) recordings made in multiwell MEA (mwMEA) plates over the first 12 days in vitro (DIV). In primary cortical cultures, action potential spiking activity developed rapidly between DIV 5 and 12. Spiking was sporadic and unorganized at early DIV, and became progressively more organized with time, with bursting parameters, synchrony, and network bursting increasing between DIV 5 and 12. We selected 12 features to describe network activity; principal components analysis using these features demonstrated segregation of data by age at both the well and plate levels. Using random forest classifiers and support vector machines, we demonstrated that four features (coefficient of variation [CV] of within-burst interspike interval, CV of interburst interval, network spike rate, and burst rate) could predict the age of each well recording with >65% accuracy. When restricting the classification to a binary decision, accuracy improved to as high as 95%. Further, we present a novel resampling approach to determine the number of wells needed for comparing different treatments. Overall, these results demonstrate that network development on mwMEA plates is similar to development in single-well MEAs. The increased throughput of mwMEAs will facilitate screening drugs, chemicals, or disease states for effects on neurodevelopment.EC was supported by a Wellcome Trust PhD studentship and NIHR Cambridge Biomedical Research Centre studentship. DH was supported by student services contract #EP-13-D-000108 and by a travelling fellowship from the Company of Biologists.This is the final version of the article. It first appeared from SAGE Publications via https://doi.org/10.1177/108705711664052

meaRtools: An R package for the analysis of neuronal networks recorded on microelectrode arrays.

Here we present an open-source R package 'meaRtools' that provides a platform for analyzing neuronal networks recorded on Microelectrode Arrays (MEAs). Cultured neuronal networks monitored with MEAs are now being widely used to characterize in vitro models of neurological disorders and to evaluate pharmaceutical compounds. meaRtools provides core algorithms for MEA spike train analysis, feature extraction, statistical analysis and plotting of multiple MEA recordings with multiple genotypes and treatments. meaRtools functionality covers novel solutions for spike train analysis, including algorithms to assess electrode cross-correlation using the spike train tiling coefficient (STTC), mutual information, synchronized bursts and entropy within cultured wells. Also integrated is a solution to account for bursts variability originating from mixed-cell neuronal cultures. The package provides a statistical platform built specifically for MEA data that can combine multiple MEA recordings and compare extracted features between different genetic models or treatments. We demonstrate the utilization of meaRtools to successfully identify epilepsy-like phenotypes in neuronal networks from Celf4 knockout mice. The package is freely available under the GPL license (GPL> = 3) and is updated frequently on the CRAN web-server repository. The package, along with full documentation can be downloaded from: https://cran.r-project.org/web/packages/meaRtools/.Wellcome Trust (+ other USA sources